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Reliable long-time storage of arbitrary quantum states is a key element for quantum information 
processing. In order to dynamically decouple a spin or quantum bit from a dephasing environment 
by non-instantaneous pulses, we introduce an optimized sequence of A'' control tt pulses which are 
realistic in the sense that they have a finite duration and a finite amplitude. We show that optimized 
dynamical decoupling is still applicable and that higher-order decoupling can be reached if shaped 
pulses are implemented. The sequence suppresses decoherence up to the order ^^(T^^^) -I- ©(r*^), 
with T the total duration of the sequence and Tmx the maximum length of the pulses. The exponent 
M G N depends on the shape of the pulse. Based on existing experiments, a concrete setup for the 
verification of the properties of the advocated sequence is proposed. 
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I. INTRODUCTION 

In quantum information processing (QIP) and in nu- 
clear magnetic resonance (NMR) it is essential to be able 
to decouple the quantum bit or the spin, respectively, 
from its environment. Both fields of research are of wide- 
spread interest and relevance. In the former the ultimate 
goal is to realize reliable long-time storage of quantum 
information with as low as possible error rates. This is a 
prerequisite for QIPl'-. In the latter, the high-precision 
measurement of nuclear spin dynamics is a long-standing 
goat^i^. 

Besides choosing well-isolated systems the application 
of appropriately tailored sequences of control pulses^i^i^, 
i.e., dynamic decoupling (DD), is one of the promising 
routes to this goal. The basic idea goes back to Hahn's 
spin echo pulse which averages a static perturbance to 
zerc^. For a dynamic environment, or bath, sequences 
of pulses are required^i^"— . The early suggestions are 
essentially periodic in time. 

Recently, the additional advantages of sequences with 
non-equidistant pulses were discovered. Concatenation 
(CDD) can suppress unwanted couplings in a high power 
T' of the length of the sequenceii. But for the method 
used in Ref. Illl the required number N of pulses grows 
exponentially with 4^. For pure dephasing, it was shown 
that this growth can be reduced decisively to a linear one 
Z oc A^ if the instants of tt pulses were chosen according 
to 

=Tsin2(j7r/(2VV + 2)), (1) 

which is called UDD (Uhrig DD). The relation ^ was 
derived for a spin-boson modeliS where it was observed 
that no details of the model entered. On the basis of nu- 
merical evidence and finite order recursion it was conjec- 
tured that ([T]) is applicable to any dephasing modeli^^^. 
This claim was finally proven^^ for any order in T. For 
various simulated classical noise spectra the experimental 
verification of the theoretical results was achievedi^""— by 
microwave control of the transition in Be ions. The sup- 



pression of the decoherence of the electron spin of hydro- 
gen radicals was investigated by electron spin resonance 
in crystals of irradiated malonic acidic. The decoherence 
was due to the quantum noise induced by nuclear spins. 
Again, the UDD proved superior to standard sequences. 

For general decoherence, concatenation of the UDD se- 
quence (CUDD) can be used2fl. For a suppression of the 
decoherence up to T^, the number of pulses grows as 2^ 
which is an improvement by a square root with respect 
to the CDD of Ref. fill. A more efficient scheme, called 
quadratic DD (QDD), which requires only a quadrati- 
cally growing number of pulses, has been proposed very 
recently based on numerical^! and analytical evidence^^. 

AU these sequences (periodic DDia, CDDii, UDDi^, 
CUDD^", and QDD^^'^S) rely on instantaneous, thus 
idealized, pulses. This problem was realized early on 
and ongoing research investigates pulses of finite dura- 
tion T^^— and sequences of such pulsesS^Ti^^. Eulerian 
P q28^33 jg (designed to annihilate the first order of a Mag- 
nus expansion over the whole sequence. Thus corrections 
of the order of Trp are not excluded. Similar caveats ap- 
ply to many other sequencea^"— . Also the experimental 
realizations in Refs. [la - llSl have to take into account that 
real pulses cannot be instantaneous because the control 
amplitudes are necessarily bounded. 

Our aim here is to derive an optimized sequence with 
UDD properties which relies on realistic pulses of finite 
duration and which is adapted to these real pulses. We 
do not provide a general scheme to use pulses of bounded 
control for arbitrary DD sequences. If the shape is ap- 
propriately designed, the pulse can be approximated as 
an instantaneous one up to 007^^). For M = 3 explicit 
results were derived in Ref. |3J while a recursive scheme 
for arbitrary M has been proposed recently'^^. As far as 
the correction ©(r*^) is negligible, the proposed sequence 
displays the same exact analytic properties as the UDD 
sequence of ideal, instantaneous pulses. 

We approach the problem hierarchically. That means 
that we use the pulses and the periods of free evolution 
as building blocks for the sequence. First, the properties 
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of the pulses are derived and discussed. Second, these 
properties are used in the sequence. To this end, we ex- 
ploit scaling in two independent variables, namely the 
durations Tp of the pulses, whose maximum is Tmx) and 
the total duration of the sequence denoted by T. Note 
that these two time scales are largely independent, both 
in theory and in experiment, because the pulses are not 
applied back to back. The only constraint is that T must 
be larger than the sum of the pulse durations Tp (be- 
longing to pulse j) 

N 

T>Y.^P.,- (2) 

Relying on pulses, which cancel all orders m < M, the 
whole sequence avoids all mixed terms T"r™^, where 
n < A'^ -|- 1 and m < M. Hence important progress over 
existing proposals^^"— is achieved. 

II. MODEL 

We start from the Hamiltonian 

H^Ao + a,AiF{t) (3) 

where az is the z component of the Pauli matrices. It is 
acting on the 5* = 1/2 spin or, generally, on the two- level 
system which represents a qubit. The operators Ai act on 
the bath only; they may also be c-numbers. We consider 
any kind of bath with bounded operators for the sake of 
the mathematical argument \ \H\ \ < 7 < 00, where || • || is 
any appropriate operator norm which remains invariant 
under unitary transforms. We expect that the order of 
suppression of the decoherence holds for any bath which 
can be approximated by bounded baths, i.e., the bath 
should have a hard high-energy cutoff. 

No spin flip terms are included in ([3]) implying an in- 
finite spin- lattice relaxation time Ti. It is an excellent 
approximation if Ti ^ T2 where T2 is the dephasing 
time. Such a situation is achieved in the rotating refer- 
ence frame of a system where the two levels with eigen- 
values ±1 of (T2 lie energetically far apart. Longitudinal 
relaxation and general decoherence will be addressed be- 
low. 

Moreover, ^ is the effective Hamiltonian in the in- 
teraction picture of the short control pulseai^rJ^. Thus 
the switching function F{t) G M appears. The simplest 
example is an instantaneous tt pulse at t = tj which real- 
izes a rotation about an axis perpendicular to az- Then 
F{t) changes its sign at t — tj abruptly while it is con- 
stant elsewhere. A sequence of such pulses at the in- 
stants {tj} with j e {1,2,..., TV} implies F{t) = (-l)J 
for t G {ij, tj+i} where we define to = and ijv+i = T. 

Our derivation is based on the bounded quantum 
model ([3]) . Thereby, classical Gaussian noise is treated to 
the extent that it can be approximated by the quantum 
model^^. Certainly, non-Gaussian classical noise, see for 



instance Ref. |H, should be considered separately which 
is beyond the scope of the present article. 

III. DERIVATION 

Optimization of the sequence means to ask the ques- 
tion which switching instants tj make the sequence 
{tj} most efficient. For ideal instantaneous pulses, it 
was shown that the UDD instants (HJ are optimum in 
the sense that the time evolution depends on the spin 
weaklyi^rJ^. The time evolution operator 

N 
3=0 

for the eigenstates of ctz with eigenvalues ±1 depends on 
the spin only in a high power of T 

U+-U^= 0((7T)^+i). (5) 

The analytical derivation of ^ is achieved by direct 
time-dependent perturbation theory (TDPT)i^ in pow- 
ers of tH. Thus, if the N + 1st power does not vanish 
it is of order (7T)^+^. The iterated time integrations 
of TDPT are conveniently expressed by the substitution 
t :— T SIT? [-0 /2) as integrations over the variable 'd. The 
instants (P) are equidistant if expressed in 'd because for 
F{d) -.^ F{Tsm'^{^/2)) we have 

with j G {0, A^}. Allowing j to take all integer val- 
ues J G Z the function i^uDD becomes an odd function 
with Fudd(i? + 7r/(iV + 1)) = -Fudd(i?)- Hence the 
Fourier series of i^uDD(^) comprises only odd sin har- 
monics sin(;(iV -f !)■&) with I G {1,3,...}. The coeffi- 
cients are 4/(7r/). From this property, ([5]) is derived by 
exploiting trigonometric addition theorems recursively^^. 

The power of the UDD sequence has been demon- 
strated experimentallyi^iii. The noise, i.e., the coupling 
to the bath Ai, is simulated, so that it can be switched 
off during the pulse. Thereby, a partial solution of the 
finiteness of the pulse amplitudes is achieved. But gener- 
ally decoherence processes cannot be switched off. Using 
pulses of duration Tp with constant amplitude instead of 
instantaneous pulses introduces an unwanted term of the 
order 7Tp at each rotation, i.e., linear in the pulse length. 
For a sequence of length N these corrections can accu- 
mulate to NjTp unless the contributions of subsequent 
pulses cancel each other. 

An improvement by one order in 7Tp is achieved by the 
ersatz n pulse which makes the linear correction vanish. 
Then the time evolution operator Up for a pulse reads 

[/p(t + Tp,t) = U;f^^\t + rp,t) + Oii-fTpf') (7a) 
U^dc^\t + Tp,t) = e-'(^p-^=)^P0e-'^=-^, (7b) 
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where M = 2. It is understood that t marks the be- 
ginning of the pulse and t + Tp its end. It is important 
that H is the Hamiltonian of the total system, i.e., spin, 
bath, and their mutual coupling, Pg is the ideal pulse 
with 9 = TT, and Tg is the instant when the approximated 
ideal pulse occurs. In a sequence {tj}, the instant Ts is to 
be identified with the switching instants tj. No adjust- 
ment of the sequence takes place. Relation ([7]) can be 
achieved by shaping the pulse appropriately^^"— . Hence 
we can set up a UDD sequence with more realistic pulses 
of the kind ([7]) for which the deviations read 



U 



UDD 



UDD 



M\ 



(8) 



The property (fTTj) is promising, but it cannot be used 
in standard DD, or in UDD in particular, as ersatz for 
an instantaneous pulse. This is so because any standard 
DD sequence presupposes that between the pulses 
the full Hamiltonian i7, not only ^o, is active. This 
conceptual obstacle cannot be solved by pulse shaping 
because the no-go theorems block further progres o^^'^^ . 
To overcome this obstacle is the main achievement of 
the present paper. We find that an adjustment of the 
sequence to the pulses of finite duration is required. 

Our present fundamental observation is that relation 
(lllbp translates to F{t) for a single realistic pulse be- 
tween t~ and = <" -|- Tp in the form 



with AI — 2. The additivity of the corrections is a 
straightforward property of the unitary evolution oper- 
ators. If we denote the UDD sequence made from the 
ideal pulses 
know from (j7ap 



ideal 



(t -I- rp,t) in dZEl) by 



UDD, ideal 



/rUDD _ AUDD, ideal 
(7± — 



0(Af(7T-mx) 



(9) 



for N pulses. The unitary invariance of the norm 7 is 
used for each pulse. Next, we know from the properties 
of the UDD sequenceH^i^ 



^yUDD, ideal _ jjUDD, ideal 



0{{lT) 



N+l\ 



(10) 



Combined with ^ this equation implies (|H]). 

The bound {'^T)^'^^ resulting from the sequence can 
be improved systematically by enlarging N . The bound 
-^(7''mx)^ resulting from the pulses can be improved by 
making it shorter. But if this is not possible, one is stuck 
because the exponent of 2 cannot be incremented for 9 = 
TT pulses as implied by mathematical no-go theorem a^^'^^ . 
Hence we are facing here a serious conceptual obstacle. 

Recently a variant of ([7]) 

C/p(i + rp,i) = U;^^°{t + T^,t) + 0{{^T^f') (11a) 

(lib) 



Uir°{t + T^,t) 



g-irpAop 



with M — was shown'' to reduce the correction to 
(7Tp)^. Note that in (fTTI) only the Hamiltonian Aq of 
the bath occurs without coupling to the spin. Hence 
[Aq, Pg] = holds and no Ts needs to be introduced. 

Explicit solutions are obtained for pure dephasing'^^. 
The correlation time of the dephasing bath should not 
be much smaller than Tp. Moreover, no no-go theorem 
was found which prevents to achieve higher orders as 
well. Indeed, a recursive scheme based on concatenation 
is proposed which achieves arbitrary order M at expo- 
nential cost^^, i.e., each composite tt pulse consists of 
> 17*^-1 elementary pulses. This demonstrates that in 
principle arbitrary M can be achieved though the expo- 
nential cost may spoil its practical usefulness. But due 
to the shortness of the pulses compared to the whole se- 
quence (Tp T) we do not expect that particularly large 
values of M are required. 




t < t- 
<t<t+ 
t > t+ 



(12) 



The TT pulse implies the inversion of the sign. But during 
the pulse itself the relation ((TT|) implies that the coupling 
between spin (qubit) and bath is effectively averaged to 
zero up to C'((7rp)*^). This is so since Aq in (jllbp does 
not comprise the spin-bath coupling; it only comprises 
the bath dynamics. This implies that the switching func- 
tion F{t) takes the value zero during the pulse. Note that 
there are jumps in the switching function even though 
the pulse is generated by bounded control. The reason 
for this behavior is that unitary time evolution is con- 
sidered over finite time intervals, not over infinitesimal 
intervals. This means that from the hierarchical level of 
the sequence we do not look into the pulses. The descrip- 
tion with F{t) is only valid on the level of the sequence, 
not within the pulse interval. Furthermore, the correc- 
tion term in Eq. (jllap may not be forgotten. 

Next we look for a sequence with F{t) G { — 1,0,1} 
which corresponds to an odd function F{{}) G {—1,0, 1} 
with the antiperiodic behavior F{'d + n/{N + 1)) — 
—F{-d). Such a sequence, RUDD (realistic UDD), allows 
for the same mathematical argument as UDD ([5]) en- 
suring that the effective time evolution of the spin is the 
identity up to corrections of the order {'-^T)^^^ . The rea- 
son is the antiperiodicity of the switching function which 
is the fundamental reason for the annihilation of the pre- 
ceding orders^'^'^^. Note that this argument holds only for 
UDD and similar optimized sequences. Hence we do not 
provide a general scheme for the incorporation of pulses 
of finite duration into arbitrary sequences. 

The sequence fulfilling the requirement of antiperiod- 
icity reads 



Fi 



RUDD 



.1). 





for d e 
otherwise 



Iv+T 



P' N+1 



(13) 

for j G The Fourier series comprises only the 

odd sin harmonics sin(^(A^ -I- l)-d) with coefficients 
Acos{l{N + l)'dp)/{TTl). The parameter < -dp < 
tt/{2N+2) determines the duration of the pulses. Except 
for the given inequality it is independent of N. Note that 
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pulses of equal duration in •& do not correspond to pulses 
of equal duration in time t 



"RUDD 



it) 



i-iy 





forte (t+,tT^O 
otherwise 



(14) 



with tf Tsin^ 



This is illustrated for 



± ,9 /2 

2N+2 ^ ^P'^ 

N — {1, 2, 3, 4} in Fig. [T] where also the necessary time- 
dependent amplitudes v{t) defining the control Hamilto- 
nian Hc(t) = v{t)ay are shown; ay is the y component 
of the Pauli matrices. For instances, the amplitudes can 
be parametrized by 

ve{t) = 9/2 + {ae - 6/2) cos(27rt/Tp) + 
[be - ae) cos(47ri/rp) + 

[c0 — bg) cos(67ri/rp) — ce cos(87rt/Tp). (15) 

There is one subtlety about the beginning and the end 
of the sequence. In order to generate the switching func- 
tion (fH|) there must be a first and a last pulse which av- 
erages the coupling between spin and bath to zero while 
inducing no net rotation. For this purpose 9 can take 
any multiple of 27r. Solving the equations derived in Ref. 
|34| . which imply that the pulse fulfills the relation ([TT|). 
leads to the parameters given in the caption of Fig. [TJ 




FIG. 1: Amplitudes v{t) (black solid and blue short-dashed 
lines) for pulses rotating about axes in the xy plane and 
the resulting switching functions F{i) (red dashed lines) for 
= 1,2, 3, 4 TT pulses. For clarity, fairly large values i^p = 
Q.7n/{2N -\- 2) are chosen for the sake of clarity. For render- 
ing purposes, i?p is chosen here to depend on N . Generally, it 
only has to fulfil i9p < n/{2N + 2). Black solid lines stand for 
9 ^ TT pulses with u,r(t) as in (fTSj) with a,r = 10.804433[l/rp], 
6^ = 6.831344[l/rp], = 2.174538[l/rp]. The first and the 
last pulse rotates by 2n (blue short-dashed lines) with a2T, 



10.236155(1/7 



= 2.9661717[l/rp], C2n = 0.889052[l/rp 



For clarity, U27r(i)/5 is plotted. 

The duration of the pulses is shortest towards the ends 
of the interval T for which the quantum state of the spin 



is to be stored. Concomitantly, the amplitudes are largest 
for the first and the last pulse. In practice, the initial and 
the final pulse can be combined with the pulses by which 
the quantum state of the spin (the qubit) is generated, 
for instance a 7r/2 pulse (for solutions see Ref. 34). 

What has been achieved by the sequence (fH)) depicted 
in Fig. [T]? This sequence is an optimized dynamic de- 
coupling scheme made from pulses of finite duration and 
finite amplitudes with analytically founded properties. 
Bounded control is a crucial aspect for realistic sequences, 
so that the proposed sequence is an important step closer 
to a realistic scenario. Nevertheless, the sequence is still 
optimized in the sense that it shares the same power law 
property as the UDD built from instantaneous pulses. 

There are two sources for corrections in the unitary 
time evolution [/^udd ^j^g RUDD sequence. The first 
kind of corrections stems from the pulses which are only 
close to Up'^'^°{t+Tp, t) but not identical to it, see Eq. pT|) . 
Denoting the time evolution of the RUDD sequence made 
from the pulses f7^"°(t + rp, t) in (ITTb)) by U. 
know from (|llap 



RUDD, zero 



we 



U. 



RUDD 



U. 



RUDD, zero 



0(iV(7r„ 



\M\ 



(16) 



for N pulses. The second kind of corrections stems from 
the sequence itself. The time evolution ^/^udd.zoio 
rigorously governed by Frudd(i^) defined in p^ . Then 
we know from Ref. \m that 



RUDD. zero 
+ 



U: 



RUDD, zero 



OihT) 



N+l\ 



(17) 



The total correction is given by the sum of both kinds of 
corrections because their norm is invariant under unitary 
transformations. So in analogy to ([8]) we obtain 



RUDD 



RUDD 



O ((7T) 



N+l\ 



(18) 

We stress that this relation excludes mixed terms 
(7T)"(7rnix)" with n < iV + 1 and m < M because 
each pulse complies with (|llap separately. We point out 
that M does not need to be as large as N because the 
pulses are much shorter anyway. So the relatively short 
and simple pulse found in Ref. |33 realizing M = 3 may 
often be completely sufficient. 

For later reference, we point out that the above deriva- 
tion also holds if we allow for an explicit analytic time 
dependence of the operators Aq and Ai in This 
was recently shown by us in the context of optimized dy- 
namic decoupling for time dependent Hamiltoniana^^ re- 
lying only on the mathematical properties of the switch- 
ing function Fij!}). Hence the same argument also applies 
to the RUDD sequence if the pulses are shaped to realize 
zero coupling during their duration, see Eq. (jl2l) . This is 
definitely the case if there is no time dependence during 
the pulses because the pulses suggested in Refs. HI and 
m can be used. This is indeed a relevant case as we will 
discuss below. 

We emphasize that the RUDD with (fT5|) provides an ef- 
ficient scheme for dynamic decoupling based on bounded 
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control. It is the main result of our paper. The pre- 
viously obstructive no-go theorem a^^'^^ can be circum- 
vented by the RUDD approach. The qualitative novel 
finding in the present work is that the sequence has to 
be adjusted in a precise way in order to allow for realis- 
tic pulses while preserving the properties of the sequence 
of ideal pulses. Above we constructed a precise prescrip- 
tion which achieves the necessary adjustment. We expect 
that this observation extends beyond the case of UDD 
and RUDD. This expectation is illustrated in the next 
section. 

We emphasize that the number of pulses N cannot be 
made infinite without using shorter and shorter pulses 
with larger and larger amplitudes. Hence a given bound 
to the available power of the control pulses limits the 
maximum possible number of pulses for a given inter- 
val T. But such limits exist in any experimental setup 
anywayiSii^iii^ and we expect that the RUDD approach 
will prove its usefulness for a moderate number of pulses. 
The limit T — is studied here to characterize the 
mathematical properties of the idealized situation. The 
achievement of the RUDD over the UDD sequence is that 
for any finite duration T and finite number of pulses N 
only pulses of finite amplitude are needed. 

IV. ITERATED SEQUENCES 

In view of the above we expect that the famous CPMG 
sequence*'^ can be improved for realistic pulses as follows. 
The CPMG is given by the n-fold iteration of the two- 
pulse cycle t — TT — 2t — TT — where tt stands for a tt 
pulse and t for free evolution of time t. This two-pulse 
cycle is the UDD sequence for A^ = 2 pulsesi^. Hence for 
pulses of finite duration the iteration of the N = 2 panel 
in Fig. [1] suggests itself. A slight modification is possible 
by replacing two 2t: pulses, where two cycles meet, by 
one 27T pulse of double the length. Hence it is promising 
to use the sequence 

(2^)t, [-t2 - TTr, - 2t2 - -Kr, - ^2 - (2^)2tJ""' - 
~t2 - vr,^ - 2^2 - -Kr, -t2- (2^)t, , (19) 

with ti = 2t{l - cos(i?p)), t2 = 2tsin((7r/6) - i9p), and 
Tp = 4t cos(7r/6) sin(i9p). The subscripts indicate the 
pulse durations. We iterate that the advocated recipe to 
account for bounded control only applies to UDD-type 
sequences. 

V. SIMULATION OF A RUDD SEQUENCE 

The advocated RUDD sequence relies on its mathe- 
matical properties which have a certain beauty in them- 
selves. But the ultimate check will be its experimental 
usefulness. A crucial step on this route is an experi- 
ment with simulated noise such as the one performed for 
UDD^^"^*. There, the simulated noise was switched off 



during the pulse. The theoretical calculations took this 
dead time of the noise into account. But variable pulse 
lengths such as in RUDD were not considered. 

We propose to implement the RUDD according to p4)) 
with pulses of finite, constant amplitudes during the in- 
tervals where F{t) = 0. No pulse shaping is required if 
the noise is switched off during the pulse so that F{t) = 
is fulfilled by construction. Hence we have 

jyRUDD _ ^RUDD ^ ^ {{jTf+') (20) 

for this particular experiment instead of (jl8p . 

The pulse intervals have to be chosen as in Con- 
comitantly the amplitudes have to vary to ensure that 
the pulses are tt pulses. In this way any deviation result- 
ing from the pulses is eliminated. It is highly interesting 
to investigate if such a RUDD sequence is more powerful 
than existing realizations. 

VI. LONGITUDINAL RELAXATION 

A UDD sequence can also suppress longitudinal 
relaxation^^. Pulses of angle tt about the z axis can sup- 
press terms proportional to x and y component, i.e., 
and ay, of the Pauli matrices up to order (-fT)^'^^ for 
{tj} as in ([T]). Concatenation of such UDD sequences 
(CUDD) can be used to suppress any kind of relaxations^. 
The QDD appears to be the most efficient scheme to ful- 
fill this purpose^ii^. 

The pulses depicted in Fig. [1] and computed in Ref. [s^ 
also work to order (7Tp)'^ if used for rotations Pg around 
the z axis for arbitrary couplings to ax, cry, and Uz ■ The 
pulse induces an inversion of the sign of the couplings 
along and ay. 

To see this one has to modify the specific calculation 
for a rotation about a fixed axis in Ref. according to 
Aq Bq — Aq + AzGz and a- A — > a^- A — axA^ + ayAy 
for a Hamiltonian H ~ Aq + a ■ A. The Eq. ([TT|) becomes 

Upit + Tp,t) = e-'^p^«F|C/G(Tp,0), (21) 

where Ug{tp, 0) encodes the corrections. It is given by 

?7g(tp,0) ^r{e-'^«^'^(*)'^*}. (22) 

where the time dependent Hamiltonian of the corrections 
G{t) stands for 

G{t) := e'^^\P^y {a^ ■ P/e-'^«* (23) 

with P^ = exp {^—i(Jz /g ds v{s)^ resulting from Hc{t) = 
v{t)az representing the pure control rotation at instant 
t. We show that (Jg{tp, 0) = 1 + 0{{-iTpf) if the shaped 
rotations about Uy proposed in Ref. |3j are applied about 
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The Magnus expansion^iSl allows us to write the time 
evolution in Eq. (1^^ in terms of cumulants Uc{Tp,0) — 
exp (— zTp ^^-^ ?7^'^) . Each cumulant Ty*-*-* scales as 
{'-fTpY. Following the approach of Ref. Is^. it is straight- 
forward to find 

v''^^ = mi i^^xAy - ayA:^ + rii2 {ctxAx + (JyAy^ , (24) 

with 7711 and 7712 the first order corrections. For the sec- 
ond order one finds 77^^' = ■q'^'^'^^ + ry^^ft) 



(2a) 



7712 := / Atcosip{t) 
'0 



ri2i[B^,(JxAy - <TyAx] 
+m2[Bo, cTxAx + cTyAy] (25a) 
r/2b) = 2m3a,{[Ay,Ax]+i{Al + Al)Y (25b) 

The expressions for 7721, 7722, and r/23 are 

rjii :— I dtamilj{t) (26a) 

(26b) 

7721 / dtt sin?/>(t) (26c) 
Jq 

7722 := I dtt cosij{t) (26d) 

Jo 

V23 ■■= jj^ dtidt2 sin(^/.(ti) - V(t2))sgn(ii - t2), (26e) 

where %l^{t) := 2 J^v{t')dt'. These conditions are exactly 
the same as those reported in Ref. [13 for pure dephas- 
ing. Hence they have the same solutions and the pulses 
depicted in Fig. 1 make the first and the second order 
corrections vanish also for longitudinal relaxation. No 
changes in the pulse shapes are required. Up to the third 
order, the transverse coupling is suppressed and only the 
z-coupling survives unaltered. Between two subsequent 
pulses the sign of the a: and y coupling is inverted. For 
pulses corrected in higher order M > 3 corrections we 
again refer to Ref. |35| for a proof-of-principle construc- 
tion. Hence, the RUDD sequence is equally applicable 
for the suppression of longitudinal relaxation. 

Eqs. ([25)1 hold generally for the suppression of deco- 
herence perpendicular to the fixed axis of rotation of the 
pulse. The decohering coupling along this axis is not sup- 
pressed. The case of pure dephasing can be seen as spe- 
cial case of the more general case discussed here: There 
is no coupling along the axis of rotation and only one 
(out of two possible) perpendicular coupling. 



axes are used. Available schemes rely on recursive con- 
catenation as for CDDii or CUDD^ or on a single step 
concatenation as for QDD2ii2^. Hence it is natural to 
consider concatenation of RUDD sequences of rotations 
about two perpendicular spin axes. 

For simplicity we consider the QDD scheme which com- 
prises two levels. On the first level two (e.g.. Ax and 
Ay) of the three couplings A to the components of a are 
eliminated up to a certain order. This is exactly what is 
achieved by a RUDD of Nz rotations about z for longi- 
tudinal relaxation as discussed in the previous section.— 
Up to the corrections O {{^Tz)^-+^) + O (A^(7r^,mx)*0 
of the primary level, the resulting time evolution is given 
an effective Hamiltonian which implies dephasing only. 
Note that is the duration of the primary RUDD se- 
quence. 

The effective Hamiltonian is of the form given in Eq. 
but with time dependent operators Ao(i) and Az{t). 
The time dependence of these operatores is analytical 
since it results from the time evolution for the time inter- 
val Tz given by the Schrodinger equation on the primary 
level. Note that it is understood that all the switching 
instants are chosen relative to Tz . Then a UDD sequence 
of duration T± can be applied on the secondary level^^ 
which consists of iV_L rotations about the spin x or y axis 
to suppress dephasing up to corrections O ((7^!)^-^+^). 
This means that general decoherence can be suppressed 
by a RUDD on the primary and a UDD on the secondary 
level. 

To obtain a quadratic scheme using pulses of finite du- 
ration we use a RUDD also on the secondary level, calling 
the resulting scheme QRUDD. This is possible since in 
the derivation of the RUDD for pure dephasing in Sect. 
IHII we mentioned that the initial Hamiltonian may dis- 
play an analytic time dependence. This effective time 
dependence results here from the pulse sequence on the 
primary level. Thus it is by construction not present dur- 
ing the secondary pulses. These secondary pulses have to 
be constructed in the presence of general decoherence a- A 
so that the explicit solution in Ref. 113 cannot be used. 
But the equations to be solved are given in sufficient gen- 
erality in this reference. For a proof-of-existence we refer 
to the work by Khodjasteh et al. where concatenated 
solutions for such pulses are constructed recursively- ^. 

Hence, the known mathematical properties of UDD se- 
quences and of TT pulses suffice to conclude that even 
general decoherence can be efficiently suppressed by dy- 
namic decoupling with bounded control by means of this 
QRUDD scheme. It is a quadratic scheme of UDD se- 
quences of bounded, and thus essentially realistic, control 
pulses. 



VII. CONCATENATION OF RUDD 
SEQUENCES 

To tackle general decoherence the combination of at 
least two sequences of rotations about perpendicular spin 



VIII. SUMMARY 

We derived in this paper that an optimized sequence 
of realistic pulses (RUDD), i.e., of finite duration and 
amplitude, can be set up which suppresses dephasing or 
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longitudinal relaxation up to T^+i in the length of the 
sequence and up to in the maximum duration of the 
pulses, avoiding all mixed terms in contrast to previous 
proposals. 

This statement is based on rigorous analytical calcu- 
lations for bounded baths and it is expected to apply 
to systems with hard high-energy cutoff. Our argument 
is based on the fundamental mathematical property of 
the optimized sequences of UDD-type, namely a certain 
antiperiodicity in the auxiliary variable "d. Thus it only 
applies to such sequences and not to arbitrary sequences. 

We introduced and exploited the concept of double 
scaling in the durations Tp of the pulses and in the dura- 
tion T of the whole sequence. We emphasize that both 
scales can be varied independently except for a certain 
constraint, see Eq. 

The key achievement is to establish a precise prescrip- 
tion how the sequence has to be adjusted to allow for the 
use of pulses with bounded amplitudes, which are thus 
decisively more realistic. Only the adjustment of the se- 
quence to the use of tailored pulses of finite duration 
allowed us to circumvent no-go theorem a^^'^^ concerning 
the properties of tailored pulses. 

The proposed RUDD can be used for suppressing pure 
dephasing, i.e., suppressing coupling of the bath to one 
spin component, or for suppressing longitudinal relax- 
ation, i.e., suppressing coupling of the bath to two spin 
components. General decoherence, i.e., suppressing cou- 
pling of the bath to all three spin components, cannot 



be suppressed by a single RUDD but by a quadratic con- 
catenated scheme (QRUDD) of two RUDDs, made from 
rotations about two perpendicular spin axes. 

Based on the known properties of UD D^^i^'^i^^ and 
QDD2ii22j we think that the design of the sequences is 
very close to its optimum. But we expect that the design 
of the pulses can still be improved. While for M = 3 
(leading non- vanishing correction is cubic in Tp) rather 
simple pulse shapes are knowr^, for higher order pulses 
with M > 3 recursive concatenation provides a recipe 
for their construction at the expense of an exponential 
increase in the number of elementary pulses^. 

Certainly, further research is called for to determine 
the performance of RUDD and QRUDD for specific mod- 
els. One important issue is to determine the size of 
the prefactors of the neglected terms. Another issue on 
the way to the experimental application of RUDD and 
QRUDD is to investigate the robustness of both the tai- 
lored pulses and the sequences to imperfections such as 
imprecise timing. 

To stimulate further research on the experimental side 
we proposed an experimental setup to verify the RUDD 
for simulated noise which can be switched ofii^"-i^ so that 
a RUDD can be checked without pulse shaping. 
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